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gas,  even  though  the  plasma  return  current  fraction  increases  rapidly  with  Ivj.  As  another 
example,  it  is  found  that  an  externally  driven  discharge  current  can  completely  eliminate 
hose  instability  in  a  fixed  conductivity  channel,  but  causes  only  a  weak  decrease  in  growth 
rate  when  the  plasma  conductivity  is  modeled  self-consistently. 
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PLASMA  CURRENT  AND  CONDUCTIVITY  EFFECTS  ON  HOSE  INSTABILITY 
1.  Introduction 

Several  recent  papers  have  developed  linear  theories  of  the  hose 
instability  of  a  pinched  charged  particle  beam,  propagating  in  a  resistive 
plasma  channel  within  a  neutral  gas.^”^  In  the  published  analytic  work  to 
date,  elaborate  models  have  been  developed  to  treat  the  beam  dynamics,  but 
the  plasma  has  been  represented  simply  as  a  fixed  ohmic  conductivity 
channel  o(r),  independent  of  time  and  axial  position  z.  For  the  case  of  a 
beam  propagating  into  initially  neutral  or  weakly  ionized  gas  the  channel 
conductivity  is  formed  by  beam  ionization  of  the  gas,  and  therefore  should 
be  treated  self-consistently  with  the  beam  dynamics.  The  significance  of 
this  linkage  has  been  understood  for  some  time,  and  two  specific  effects 
have  been  pointed  out  by  Briggs,  Lee,  and  co-workers^  for  the  case  of  beam 
propagation  into  neutral  gas.  First,  the  equilibrium  conductivity 
oo(r,5)  increases  monotonically  with  distance  ?  5  v^t  -  z  behind  the  beam 
head,  in  the  forward-  part  of  the  beam  where  recombination  and/or  plasma 
cooling  have  not  yet  established  a  balance  with  beam-driven  ionization  and 
heating.  It  Is  thus  Incorrect  there  to  treat  C  as  an  ignorable  coordinate 
and  Fourier  analyze  the  mode  dependence  on  q.  Since  the  growth  rate  of  any 
resistive  instability  is  largest  where  is  small,  the  effects  due  to 
increasing  oQ  are  Important.  The  second  effect  is  seen  when  there  is  an 
appreciable  equilibrium  plasma  return  current  I  .  If  I.  flows  in  a  fixed 

p  p  - 

conductivity  channel,  an  absolute  Instability  results  from  the  magnetic 

repulsion  between  1^  and  the  currant  1^  of  a  transversely  displaced  beam. 

However  if  the  conductivity  channel  is  formed  by  the  beam  Itself,  the  head  of 

the  beam  cannot  be  displaced  from  the  channel,  which  converts  the  instability 

into  a  convective  one  (provided  the  effect  of  phase-mix  damping  among  particles 

of  different  betatron  frequency  is  modeled  correctly1*^).  Furthermore,  the 
Mwiesertpl  approwd  October  11, 1963. 


channel,  and  thus  I  as  well,  tend  to  follow  the  hose  distortion  of  the  beam, 

P 

thereby  further  reducing  the  destabilizing  effect  of  spatial  separation 

between  1^  and  1^.  The  consequences  of  these  effects  have  been  worked  out, 

in  the  low  frequency  limit  and  elsewhere  in  terms  of  certain  parametric 

representations,  by  the  Livermore  group. ^  However  hose  instability  theory 

has  never  been  treated  in  the  straightforward  form  of  a  dispersion  relation 

analysis,  valid  for  all  frequencies,  that  includes  self-consistent 

treatment  of  the  channel  conductivity  and  the  plasma  current. 

In  the  present  paper  we  include  beam-impact  ionization  and  plasma 

recombination,  treated  self-consistently ,  in  our  conductivity  model.  These 

are  usually  the  dominant  processes  for  the  case  of  a  weakly  ionized  gas, 

except  for  a  short  region  near  the  "pinch  point"  at  the  front  of  the  beam, 

7  ft 

where  an  inductive  S  spike  can  cause  strong  avalanche  breakdown.  We 

z 

further  simplify  the  conductivity  model  by  assuming  the  electron  mobility 
is  constant,  i.e.,  independent  of  variations  in  plasma  electron 
temperature.  The  conductivity  then  depends  only  on  the  electron  density. 
Although  this  is  a  reasonably  accurate  approximation  in  many  cases,  we  have 
recently  discovered  that  a  number  of  Interesting  beam  propagation  effects 
can  be  associated  with  the  temperature  dependence.  These  will  be  discussed 
in  a  future  publication. 

We  also  specialize  to  the  limits  of  an  idealized  "beam  body”, 
discussed  in  Sec.  II,  where  1^(0  has  reached  a  plateau,  avalanche  is 
assumed  to  be  unimportant,  and  the  plasma  electron  density  n@  is  not  yet 
large  enough  for  recombination  to  be  important,  and  a  "beam  tail", 
discussed  in  Sec.  Ill,  where  ionization  and  recombination  have  reached  a 
balance  for  the  beam  equilibrium  (but  not  necessarily  for  the 
perturbation).  We  then  show  that  the  Inclusion  of  self-consistent 


$1 


M 

I 

VI 


conductivity  physics  modifies  the  problem  in  a  way  that  is  mathematically 
simple  but  dramatically  changes  the  physical  results.  Given  any  beam 
dynamics  model,  and  the  mathematical  machinery  to  solve  for  a  dispersion 
relation  in  the  fixed  channel  case,  there  is  no  further  difficulty  in 
solving  for  a  complete  dispersion  relation  with  a  self-consistent  treatment 
of  the  channel.  In  fact,  a  dispersion  relation  is  found  in  closed  form 
when  the  relatively  simple  "spread-mass”  model^*^  of  beam  dynamics  is  used. 

For  relativistic  electron  beams,  it  is  convenient  to  use  z  and  5  as 
independent  variables,  rather  than  z  and  t.  In  the  beam  body,  normal  modes 
are  found  with  5  and  z  dependence 


exp[-i<s*n(c/c  )-ifl2/v  ] 
o  z 


3  (C/CQ)  1  {cos[a»rin(c/5o)-i  sin[urin(c/i;o)  ]}exp(-iflz/vz) 


rather  than  the  usual  form 


exp(-io»c/v  -  iJJz/v  ) 
z  z 


for  a  fixed  channel.  In  Eq.  (la),  is  an  arbitrary  index  point.  In  the 
beam  tail,  the  form  (lb)  holds,  since  5  is  an  ignorable  coordinate  in  the 

equilibrium. 

Plasma  current  effects  fall  into  three  different  regimes.  (1)  If  there  is 
no  significant  equilibrium  plasma  current  (lip I  «  I^)»  conductivity 
perturbations  play  no  role.  In  the  tall,  u(n)  is  identical  to  the  dispersion 
relation  for  a  fixed  channel.  In  the  body  u(ft)  is  proportional  to  the 
frequency  »(Q)  that  would  be  found  in  a  fixed  channel.  In  fact,  for  any 


conductivity  profile  of  the  fora  oq(t,z)  •  <j  '(r)o"'(;),  the  exact  r 

,  r 

dependence  of  the  normal  modes  is  exp  -i ;’d5''j[Q;a  (r,;)],  where  Eq.  (lb)  with 
the  dispersion  relation  (r)  ]  is  the  normal  mode  in  a  fixed  channel  with 

conductivity  profile  oQ(r).  (2)  If  there  is  a  high  degree  of  current 
neutralization  (0  <  1^  +  I  «  1^),  I  is  strongly  destabilizing  but 
conductivity  perturbations,  which  allow  the  beam  to  follow  the  channel, 
partially  cancel  the  destabilizing  effect.  (3)  If  I  flows  parallel  to  1^ ,  as 
it  could  for  example  if  an  external  discharge  is  applied  to  the  channel,  it 
is  possible  in  a  fixed  channel  to  stabilize  the  hose  mode  completely,  but 
conductivity  perturbations  restore  the  instability  by  causing  the  discharge 
current  to  follow  the  unstable  beam,  rather  than  guiding  it  along  an 
established  path.  In  all  cases,  the  perturbed  conductivity  tends  to  push 
the  hose  growth  rate  toward  the  value  it  would  have  in  the  absence  of 
plasma  current.  In  particular,  we  find  that  increasing  the  beam  current, 
thereby  increasing  both  the  fractional  current  neutralization 
(destabilizing)  and  the  conductivity  a  (r,5)  (stabilizing),  usually  leads 
overall  to  a  reduced  instability  growth  rate  u>^  in  the  beam  body. 

The  tail  of  a  beam  propagating  with  a  large  return  current  in  a  weakly 
ionized  plasma  channel  is  also  subject  to  another,  even  more  important 
stabilizing  effect:  the  balance  between  impact  ionization  and 
recombination  leads  to  a  channel  broader  than  the  beam,  so  that  most  of  the 
return  current  flows  outside  the  beam  and  has  little  destabilizing 
effect.  As  a  result  of  this  and  the  effects  discussed  in  the  previous 
paragraph,  we  find  that  increasing  the  beam  current  strongly  reduces  the 
growth  rate  in  the  tail. 

In  the  present  paper,  we  do  not  include  other  sources  and  sinks  in  the 
conductivity  model,  such  as  avalanche  breakdown,  beam  heating  of  the 


plasma,  or  plasma  cooling  due  to  radiation,  conduction,  etc.  These  effects 
could  be  included  formally,  but  would  prevent  solution  for  the  ;  dependence 
in  the  closed  forms  (la)  or  (lb).  They  are  therefore  best  treated  within  a 
numerical  scheme  that  solves  for  the  r,  dependence,  e.g.  the  codes  VIPER^ 
(multi-component  beam  dynamics  ),  SIMM1  (linearized  beam  particle 
simulation),  VALIUM3  (linearized  Vlasov  solver  for  beam  dynamics),  or 
EMPULSE^  (spread  mass  model  beam  dynamics).  In  a  fully  ionized  channel, 
the  effects  of  plasma  heating  and  cooling  on  hose  instability  can  be 
expected  to  be  similar  to  those  of  ionization  and  recombination  in  a  weakly 
ionized  channel. 

In  a  previous  paper  we  have  discussed  the  effects  of  self-consistent 
conductivity  evolution  on  the  sausage  mode.  The  effect  there  was 
generally  similar  but  much  stronger;  channel  perturbations,  by  neutralizing 
the  destabilizing  effect  of  plasma  return  current,  completely  stabilized 
the  mode  over  a  wide  parameter  range. 


5 


Instability  of  the  Beam  Body 
A.  Assumptions  and  Formalism 


We  assume  that  the  beam  is  highly  relativistic  (y  >>  O,  paraxial 


(v^  «  v^  for  all  beam  particles),  and  therefore  that  vz  =>  c  for  all  beam 


particles.  The  equilibrium  beam  profile  J,  (r,-)  is  treated  as  time- 


independent.  In  the  beam's  body  and  tail,  we  assume  1^  has  reached  a 


constant  value.  We  may  also  assume  that  the  beam  radius  a(;),  defined  by 


Xb  *  *a  Jbo(t  ’  °>’ 


and  the  plasma  current  1^  are  independent  of  5.  These  latter  assumptions 


are  reasonable  because  resistive  decay  of  the  plasma  current  with  c  occurs 


slowly,  on  the  monopole  decay  scale ' 


CTo(0 


2vo  (r  -  0,  C)a  ,b< 

- e - 


where  b  is  a  large  radius  characterizing  the  region  of  space  charge 


neutralization  (b  »  a),  while  instability  growth  is  characterized  by  the 


1-3  13 

much  shorter  dipole  decay  5  scale,  * 


iroo(0«?)a 


7  8 

Furthermore,  the  beam  body  is  defined  to  lie  behind  the  pinch  point  *  ,  so 


the  conductivity  in  the  body  and  tail  exceeds  the  value  <3  *  c/4wa  required 


for  space  charge  neutrality,  and  with  the  use  of  the  paraxial  approximation 


Max  *”'s  juations  reduce  to  Ampere's  law  for  the  axial  component  A  of 


vector  potential1*, 


1 


7,2A  ,  il  J  iZL£|A. 

i  c  2  c  b  c  3? 


We  define  the  beam  body  as  the  region  behind  the  pinch  point,  but 

where  recombination  (which  depends  quadratically  on  n  )  is  not  yet 

Important.  Furthermore,  we  neglect  avalanche  ionization  in  the  beam  body; 

this  is  usually  reasonable  although  avalanche  may  well  be  important  in  the 

vicinity  of  the  E  spike  at  the  pinch  point,  ahead  of  the  body.  In  a 
z 

weakly  ionized  gas,  the  electron  mobility  is  determined  by  electron-neutral 

collisions  which  depend  weakly  on  temperature;  since  the  dependence 

of  <j  on  ng  is  much  more  important,  we  neglect  the  temperature  dependence. 

In  addition,  we  assume  that  beam-impact  ionization  is  local,  i.e.  we 
neglect  any  spatial  spreading  due  to  high-energy  secondary  electrons.  We 
may  then  write  a  simple  equation  for  the  conductivity. 


g 

jjr  o(r,c,z)  -  «Jb(r,?,z), 


where  is  the  beam  current  density  and  <  is  here  assumed  to  be  a  constant 
for  any  gas,  proportional  to  the  beam-impact  ionization  frequency  and 
Inversely  proportional  to  the  plasma  electron-neutral  collision 
frequency.  In  equilibrium,  the  solution  of  Eq.  (6)  is 


o  (r,c)  -  kJ  (r)(c-c  )  +  o  (r,c  ), 

0  DO  0  o  o 


where  is  the  front  of  the  beam  body.  The  first  term  of  Eq.  (7)  leads  to 
a  conductivity  profile  of  the  same  form  as  the  beam  profile.  We  assume 
that  the  second  term  either  is  negligible  or  has  the  same  profile,  and 
absorb  it  by  choosing  c  so  that 


A  pre-formed  channel  narrower  or  wider  chan  the  beam  cannot  be  treated 


within  the  simple  analytic  calculation  performed  in  this  section. 


Qualitatively,  a  narrower  channel  decreases  the  decay  lengths  (3)  and  (4) 


and  accentuates  the  effect  of  I  on  stability^,  which  is  stabilizing 


if  I^/I^  >  0  and  destabilizing  if  I^/I^  <  0.  The  opposite  is  true  for  a 

wider  channel. 

Since  z  is  an  ignorable  coordinate  of  the  equilibrium,  we  may  look, 
for  normal  modes  with  z-dependence  exp(-i~z/c).  We  also  follow  the  usual 


procedure  of  expanding  the  azimuthal  dependence  of  A  and  J  in  a  Fourier 


imd 


series  e  ,  keeping  only  the  monopole  (m  «  0)  and  dipole  (m  *  1)  terms, 
and  linearizing  the  dipole  terms  which  are  treated  as  small 


perturbations.  All  perturbed  quantities  ’y ( r  ,  z,  ;)  thus  take  the  form 

'•l 


z,c)  -  exp(-iftz/c  +  i9). 


Equation  (5)  for  the  dipole  A^  becomes 


(LILr. 

^3r  r  9r 


bo  r  i_U 

c  ?  9  5  i 


-  —  J 


4tt 


bl 


9A 

c 

J1  9; 


(9) 


The  second  term  on  the  right  hand  side  of  Eq.  (9)  represents  the 
current  driven  in  the  perturbed  channel  by  the  unperturbed  electric 
field  ■  -  9A  /3t.  E--  can  be  expressed  in  terms  of  I  and  a  , 


(10b) 


PWCTPP  CT  Kf  •’'7  r  •■*  <•  ^  '■'  V  V  '• 


■  .v.' 


i.e. 


J  (r)  I 

zou;  <J.  (r)5  <I.c* 

b  °  7  8 

where  we  have  neglected  the  weak  logarithmic  dependence  ’  of  E^o  on  r,  so 
that  OpQ(r)  has  the  same  radial  profile  as  0^Q(r)  and  oq( r).  Using  Eq. 
(10b),  Eq.  (9)  becomes 


t o  1  d  DO  ,  d  1, 

l3?  7  37  r  c -  ’  37^1 


4tt  j  _  4tt_  Tp  ql 

c  bl  c  <1  ; 


(ID 


Next  we  use  the  linearized  version  of  Eq.  (6),  in  the  form 


3cri  a  ai 

i  =  (l  +  5  )  —  -  icJ,  ,  , 

3?  L  s  3?'  5  bl’ 


(12) 


to  eliminate  o^  from  Eq.  (11): 

13  .  4l"tJbo<r)  .  3 

o  +  4  35)  ^3r  r  3r  r  c  5  3c^Al 


r-  (1  +  ‘*c+t?Ki- 


(13) 


Since  Eq.  (13)  has  constant  coefficients  with  respect  to  the 
operator  5(3/3?) »  solutions  exist  of  the  form 


A1(r,5)  -  A(r)  exp[ -iuj£n(? /?Q)] 


~  (jj. 

-  A(r)  ?  {cos[u>rin(?/?o)j-  i  sin[u>r£n(?/;o)]}  •  (14) 


Equation  (13)  thus  reduces  to 


9 


,d  1  d  ,  i’i«Jbo(r)>  . 

dr  f  dr  r  + -  A(r) 

c 


—  1  +  -P  V,  J.  (r).  (15) 


When  dimensionless  variables 


r  =  r/a. 


Jbo'c)  ■  •,bo(cWbo<°> ' 


are  used,  Eq.  (15)  can  be  rewritten  in  the  fora 


(^:Tr'  +  8i]jAJbo(r)  }A(r) 
dr  r  dr 


-|l  +-L_]}  (;)a2>  (18) 

1  -  iui 


which  depends  only  on  the  parameters 


f  5  Xp/Xb 


dTi  <Ib 

x  =  c  ir 


By  way  of  comparison,  Ampere's  law  takes  the  fora 


[i.  I  i_  ?  +  Siu^l  (r)]A(r)  -  -  Jfe(?)a2 
dr  r  dr  '  c 


wsfosrsst 


s  f  * 
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in  a  fixed  channel,  where  perturbed  quantities  are  of  the  fora 

A^r,;)  -  A(r)e_iuJ?/c.  (22) 

-  O 

Equations  (20)  or  (21)  must  be  supplemented  by  an  expression  for  J  (r)a“  as 

D 

a  function  of  2  and  f  and  a  functional  of  A(r);  as  we  shall  see,  this 
calculation  is  the  same  for  both  cases.  Thus  the  hose  eigenvalue  problem 
in  a  self-consistent  conductivity  channel  differs  from  the  problem  in  a 
fixed  channel  in  only  two  ways:  first,  eigenfunctions  are  of  the  form  (14) 
rather  than  (22),  due  to  the  linear  increase  of  the  unperturbed 
conductivity  aQ  with  z,  in  the  beam  body,  and  second,  there  is  an  extra  term 
on  the  right  hand  side  of  Eq.  (18),  due  to  the  inclusion  of  perturbed 
conductivity.  Since  several  successful  methods  have  been  found  to  treat 
the  problem  in  a  fixed  channel,  extension  of  these  calculations  to  Eq.  (18) 
is  immediately  possible. 


3.  Calculation  of  the  Perturbed  3eam  Current 

Since  all  particles  have  the  same  axial  velocity  =  c,  the 
coordinate  ;  is  a  constant  of  the  notion  and  plays  no  role  in  the 

A 

determination  of  J^(r),  i.e.  the  form  of  the  ;  dependence,  Eq.  (14)  or 
(22),  and  the  value  of  a  or  u>  do  not  occur  in  the  particle  dynamics. 

A 

J^(r)  does  depend  on  the  time  history  of  the  driving  potential 

mm  Q2  /  C 

A(r)e  4  (z  plays  the  role  of  time  in  the  beam  frame),  and  also  on 

equilibrium  quantities  such  as  f  and  the  on-axis  betatron  frequency 

1/2 

ngQ  =  [2T(l+f)eJ^o(0)/(ymc) J  ,  where  e  and  m  are  the  electron  charge  and 

A 

mass,  and  y  is  the  relativistic  factor.  r )  may  be  calculated  by  solving 

the  linearized  Vlasov  equation^,  by  particle  simulation1®,  or  by  using 

■> 

either  of  two  simplified  beam  dynamics  models,  the  multi-component  model 
or  the  spread  mass  model1.  The  latter  yields  a  closed-form  dispersion 
relation  which  is  useful  for  demonstrating  general  properties,  and  which 
the  reader  can  easily  solve  for  any  case  of  interest  to  him.  The  other 
dynamics  models  are  more  accurate,  but  dispersion  relations  are  obtained 
only  through  elaborate  numerical  solutions  of  the  eigenvalue  problem.  We 
have  developed  codes  that  generate  these  solutions,  and  will  cite  some 
examples. 


1.  Spread  Mass  Model 


I 

I 

> 

' 

i 


i 


The  spread  mass  model  of  E.  ?.  Lee^*^  is  the  simplest  beam 
dynamics  model  that  includes  the  essential  feature  of  phase-mix  damping  due 
to  betatron  frequency  spread.  According  to  the  spread  mass  model  for  the 
case  of  a  beam  with  3ennett  profile, 


(r) 


aQ(r,0 

0,(0, 0 


1 


?  i  o  * 

(1  +  rVaV 


(23) 


J.  ( r , ; )  and  A(r,c)  are  represented  as  rigid  sideways  displacements  of  the 

D 

equilibrium  profiles 


Jb(r.O  -  -  Y(0(dJb<j/dr)  (24a) 

A(r,c)  -  -  D(;)(dAo/dr),  (24b) 

1  14 

and  Y(c)  is  given  as  a  function  of  D(c)  by  9 

Y(c)  -  D(?)[l  +  G(Si2/nSo2)]  (25) 

where 

G(x)  =  6x{-j  -  x  +  (x-x^)[iri+£n(-i~.)  ] }  for  |x|  <  1.  (26) 

The  function  G(x)  is  plotted  in  Fig.  1  for  real  x.  Since  the  r  dependence 
of  all  quantities  is  specified  in  closed  form,  Eqs.  (13),  (23)-(26) 
loMdiately  reduce  to  a  dispersion  relation 
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( 1-iooA)  (  1 


— r1  a+f)  - 1  +  g;  ■ 


l-iui 


2  * 


'30 


(27) 


which  yields  the  solution  for  ui(2) 


U)A  ■  ”  TC  + 


f  - 


G<a 2/n  2)  ,  f-o(r.2/n  2)  2  ,  , 

- -  -)  t  T  !(>  +  t  ^  f  3  )  +4xG(a2/a..2)’l/2 


1  +  f 


3o 


(28) 


We  note  that  for  each  value  of  Q  there  are  two  solutions  for  jj,  but  one  of 
these  branches  is  found  to  be  damped  for  all  real  values  of  2,  and  thus  of 
no  interest.  Of  course  Eq.  (27)  may  also  be  solved  numerically  for  .'I(uj). 

We  may  compare  Eq.  (28)  with  the  well-known  dispersion  relation  for  a  fixed 
channel1*1^. 


2.  Multi  Component  Model 


3 

The  multi-component  model  of  beam  dynamics  calculates  the  r 

A 

dependence  of  J^,  rather  than  assuming  a  form  such  as  Eq.  (24),  and  it 
Includes  radially-localized  resonances  and  other  important  properties  of 
the  exact  solution.  According  to  Eqs.  (57),  (58)  and  (11)  of  Ref.  3,  the 
model  yields  an  expression 


Jb(r) 


2r 


1  +  f(0) 


JR  ,d  1 

Jr  R2  ^dR  R  dR  ; i  “  2 

“So 


2  -0 


3o 


dr"r'“  A(r'), 
(30) 


where 


nc2(R) 


-7T-  fo  dr  ”  tt)  \o(r)  fr^rioj1* 

K  K 


(31a) 


and,  for  generality,  we  have  defined 


t(r)  =  ■Jpo<r>/Jbo<r)- 


(31b) 


Since  J.  (r)  and  J  (r)  have  the  same  radial  profile  in  the 

DO  pO 

beam  body,f(r)  -  f  and  we  may  rewrite  (31a)  as 

,  4SJfln  o  _2  _ 

°c  <R)  5  -r-  /o dr  r  C1  -  *?)  Jbo<r)-  (3ic) 

R  R 

Using  Eqs.  (18),  (20)  and  (30),  the  eigenvalue  problem  may  be  expressed,  in 
analogy  with  Eq.  (61)  of  Ref.  3,  as  a  fourth-order  ordinary  differential 
equation 


1  (d  1-.4  -3  _  ,  d  bo  d 

t  l-r  rj  r  r  +  lux  -—  - 

8  dr  r  dr  rJ  dr 


r2  r 
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(32) 


1  Id 


1  "bo  V(0) 


with 


-  1  +  f  -  -  -  ->  i 

1  -  ij  r  dr  r  dr  2,"(r)  -  3“ 


T(r)  = 


dr'  r'“  A(r') 


and  boundary  conditions 

r(0)  -  r(0)  -  r'(-)  -  r"(-)  -  o.  (34; 

This  differs  from  the  fixed-channel  case  only  through  the  extra  factor 
[i  +  f(l  -  iu>)  M  on  the  right  hand  side  of  (32).  Thus  a  code  which  we 
developed  to  solve  this  complicated  eigenvalue  problem  in  the  fixed- 
channel  case  can  also  be  used  to  find  'ju(ft)  in  the  present  case. 


*7>;v  *>  •  y;.-  v>. 


**>■  ** 


3. 


Vlasov  Calculation  of  J. (r) 

D 


If  the  beam  equilibrium  is  taken  to  be  tiae-independent,  and 
weak  scattering  of  the  beam  off  the  gas  is  neglected,  an  exact  expression 

A  A 

for  J^r)  in  terms  of  A(r)  can  be  obtained  by  integration  of  the  linearized 
Vlasov  equation  along  unperturbed  beam  electron  orbits.  This  procedure 
yields 


3„<r)  -  -e2  /  d3e  <*(t)  +  t  dz'  -  13(e**')-’>. 


(35) 

where  p  is  the  particle  relativistic  momentum,  p  and  p,  are  the  axial  and 
perpendicular  components  of  g,  f^P^Pj^  Is  the  beam  distribution 
function,  and  (r,0,z)  are  the  particle  cylindrical  coordinates. 

9"(z'')  and  r '(z')  are  the  unperturbed  particle  trajectories  backward  in 
time,  l.e.  9"  and  r"  are  the  coordinates  a  particle  had  when  it  was 
at  z't  if  that  particle  has  coordinates  6  and  r  when  it  is  at  z. 

Thus  Q'(t')  and  r-(z')  are  also  functions  of  g  and  of  the  unperturbed 
potential  Aq. 

Equations  (35)  and  either  (18)  or  (21)  together  constitute  a  very 
complicated  integral-differential  eigenvalue  problem.  A  code,  known  as 
VALIUM,  has  been  developed^  with  the  capability  to  solve  the  fixed-channel 
case  (35)  and  (21)  as  well  as  the  self-consistent  channel  case,  Eq.  (35) 
and  (18). 
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C.  Dispersion  Relations 


In  this  section  we  discuss  the  growth  spectra  that  result 

from  the  dispersion  relation  derived  in  closed  form  in  Sec.  IIA,  or  from 
solution  of  the  eigenvalue  problems  framed  formally  in  Secs.  113  and  IIC. 
The  discussion  will  be  organized  in  terms  of  the  value  of  f  -  either  no 
plasma  current  (f»0),  extensive  current  neutralization  (0  <  1  +  f  <<  1) ,  or 
plasma  current  flowing  parallel  to  the  beam  current  (f  >  0),  and  in  terms 
of  the  effect  of  beam-generated  perturbations  to  the  conductivity,  which  is 
small  if  X  <  1  (low  beam  current)  and  large  if  X  »  1  (high  beam 
current).  The  value  of  X  depends  only  on  1^,  but  the  value  of  f  can  depend 
on  the  properties  of  the  beam  head  (rise  time  of  I,  and  avalanche 

D 

ionization  near  the  pinch  point),  and  the  presence  of  a  pre-formed 
conductivity  channel  or  an  externally  imposed  Ez  field,  as  well  as  on  Ib» 

We  conclude  by  surveying  the  case  of  beam  injection  into  field-free  neutral 
gas,  under  conditions  where  both  X  and  f  depend  primarily  on  the  single 
parameter  L^. 

We  ahall  use  the  spread  mass  model  dispersion  relation,  Eq.  (28), 
to  derive  various  analytic  results  in  closed  form.  We  also  show  plots  of 
the  dispersion  relation  for  a  variety  of  cases,  Figs.  2-5,  7  as  calculated 
from  the  multi-component  model,  Eqs.  (32)— (34 ) ,  which  is  more  accurate  and 
which  is  not  immediately  accessible  to  the  reader.  The  numerical  values  of 
the  key  features  of  the  dispersion  relations  -  maximum  growth  rates  and 
cut-offs  -  are  given  in  Table  1  for  each  of  these  models  and  for  the  Vlasov 
calculation,  which  is  theoretically  exact,  but  subject  to  numerical 
uncertainties  of  perhaps  a  few  percent.  We  find  that  both  models  agree 
well  with  the  exact  result  as  regards  the  peak  growth  rate  Max  The 
multicomponent  model  also  gives  an  accurate  estimate  of  the  width  of  the 


V 
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Table  1 

Comparison  of  the  key  features  of  the  hose  dispersion  relation  in  the  bean 
body, as  given  by  the  spread  mass  (SM)  model,  the  multicomponent  (MC)  model, 
and  a  Vlasov  (V)  calculation.  Here  Max  is  the  peak  power-law  growth 
exponent  for  growth  in  ;  and  X  Max  u  is  a  natural  scaled  quantity  from  the 
theory,  is  the  value  of  Q  for  which  *  Max  Qc  is  the  largest 

value  of  for  which  >  0,  and  Q0q  is  the  on-axis  betatron 

frequency.  These  values,  for  the  multicomponent  model,  can  also  be  read 
directly  from  the  dispersion  relation  curves  of  Figs.  2-5. 


0.044 


f 

X 

Max 

a 

m  30 

c  30 

SM 

MC 

V 

SM 

MC 

V 

SM 

MC 

V 

1.0 

-0.17 

-0.17 

0.52 

0.53 

0.5 

0.12 

0.11 

0.52 

0.55 

0.62 

0.66 

0 

0.69 

0.66 

0.52 

0.61 

0.71 

0.78 

-0.25 

1.25 

1.23 

0.52 

0.62 

0.73 

0.82 

-0.5 

2.4 

2.4 

0.52 

0.63 

0.76 

0.87 

-0.7 

4.6 

4.9 

0.52 

0.67 

0.78 

0.90 

-0.9 

15.9 

18.7 

0.52 

0.79 

0.79 

0.96 

1.0 

0.04 

0.04 

0.38 

0.42 

0.46 

0.46 

0.5 

0.16 

0.14 

0.51 

0.53 

0.63 

0.67 

0 

0.69 

0.66 

0.74 

0.52 

0.61 

0.6 

0.71 

0.78 

0.8 

-0.25 

1.22 

1.21 

1.34 

0.52 

0.62 

0.65 

0.73 

0.82 

0.9 

-0.5 

2.3 

2.4 

2.5 

0.52 

0.63 

0.66 

0.75 

0.87 

0.9 

-0.7 

4.6 

4.8 

5.0 

0.52 

0.67 

0.68 

0.78 

0.90 

0.9 

-0.9 

15.8 

18.7 

17.2 

0.52 

0.79 

0.7 

0.79 

0.96 

0.9 

1.0 

0.32 

0.29 

0.51 

0.51 

0.68 

0.72 

0.5 

0.40 

0.36 

0.42 

0.51 

0.53 

0.51 

0.69 

0.73 

0 

0.69 

0.66 

0.74 

0.52 

0.61 

0.6 

0.71 

0.78 

0.8 

-0.25 

1.12 

1.11 

1.21 

0.52 

0.62 

0.61 

0.73 

0.82 

-0.5 

2.3 

2.2 

2.5 

0.52 

0.63 

0.62 

0.75 

0.86 

0.9 

-0.7 

4.3 

4.5 

4.7 

0.52 

0.67 

0.63 

0.77 

0.89 

-0.9 

15.5 

18.3 

18.9 

0.52 

0.80 

0.71 

0.78 

0.95 

1.0 

0.64 

0.60 

0.54 

0.62 

0.75 

0.82 

0.5 

0.65 

0.62 

0.53 

0.62 

0.74 

0.80 

0 

0.69 

0.66 

0.74 

0.52 

0.61 

0.6 

0.71 

0.78 

0.8 

-0.25 

0.74 

0.70 

0.82 

0.52 

0.55 

0.56 

0.68 

0.77 

-0.5 

0.88 

0.87 

1.11 

0.52 

0.55 

0.58 

0.65 

0.78 

0.9 

-0.7 

1.7 

2.1 

2.3 

0.52 

0.67 

0.65 

0.67 

0.85 

-0.9 

11.9 

14.2 

14.8 

0.52 

0.78 

0.72 

0.76 

0.95 

0.9 
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unstable  spectrum  and  the  value  2^  for  peak,  growth  [jj^(2  )  *  Max  ^  ] ,  but 
the  spread  mass  model  underestimates  both  the  spectral  width  and  the  value 


I  -  0 
_P _ 

In  the  case  of  no  equilibrium  plasma  current,  I  ■  f  ■  0,  Eq.  (13) 
for  Xa(3)  in  the  beam  body  becomes  formally  identical  to  Eq.  (21),  which 
determines  T.a(3)  in  a  fixed  conductivity  channel.  Thus  the  dispersion 
relation  for  power-law  growth  in  the  beam  body,  Xa(3),  is  identical  to  the 
dispersion  relation  for  exponential  growth  in  a  fixed  channel  with  f  «  0. 
Perturbed  conductivity  plays  no  role  when  f  *  0,  since  the  absence  of 
equilibrium  plasma  current  indicates  that  Ezo  *  0  and  thus  a  dipole 
conductivity  a,  yields  no  dipole  plasma  current  E  a,  .  '.Then  the  spread 

*  ZO  1 

mass  model  is  used,  the  dispersion  relation  is  given  by  Eq.  (29), 


Xa(3)  -  iG(a2/3So2),  (36) 

and  is  shown  in  Fig.  1.  For  the  spread  mass  model,  the  peak  growth  rate  is 

Max  ^  -  0.69  x  1  at  3/3gQ*  0.52,  and  the  unstable  spectrum  runs  from  3*0 

to  3  *  0.71  3^q.  The  equivalent  result  from  the  multicomponent  model  is 

shown  as  the  f  *  0  curve  in  Fig.  2;  the  peak  growth  rate  is 

Max  u)^  ■  0.65  X  ^  at  3/3  q  ■  0.61,  and  the  unstable  spectrum  runs 

from  3  ■  0  to  3  ■  0.78  3„n.  We  note  that  a. (0)  ■  0,  an  exact  result. 

g0  i 

Low  Beam  Current  (  X  *  0) 


Next  we  consider  the  limit  X  ♦  0,  which  corresponds  to  low  beam 
current,  as  seen  in  Eq.  (20).  For  example,  in  air  we  use  the  estimate^ 


where  1^  is  in  kiloamperes .  Equation  (37)  corresponds  to  <  =  3.3  x  10-i4 
cm/statcoul  in  Eq.  (6). 

In  the  limit  X  -*■  0,  for  any  given  value  of  f  and  for  any  beam  dynamics 
model  the  dispersion  relation  for  \u(2)  is  identical  to  the  dispersion 
relation  for  in  a  fixed  channel.  This  is  because  the  conductivity 

changes  slowly  in  this  limit,  and  thus  the  effect  of  conductivity 
perturbations  is  negligible.  For  example,  the  spread  mass  dispersion 
relation,  Eq.  (28),  reduces  to^*^ 

f  -  G(22/Q  2) 

Xu  «  -  i  - . £  •  (33) 

If  f  <  0,  there  is  an  unstable  mode  at  a  -  0  with  growth  rate 

^(0)  »  X  l[-f(l  +  f)  *"],  an  exact  result,  but  the  most  unstable  mode  is 

still  at  n/Qa  *  0.52,  where  co  ■  Max  u>,  *  X  1  (-f  +  0.69)(1  +  f)  1.  The 
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equivalent  result  from  the  multi-component  model  is  shown  in  Fig.  2  for 
several  values  of  f.  The  peak  growth  rate  for  any  given  value  of  f  is 
close  to  that  given  by  the  spread  mass  model,  as  seen  in  Table  1,  but  the 
value  of  0  for  peak  growth  increases  with  the  return  current  fraction,  and 
the  unstable  spectrum  given  by  the  multi-component  model  is  considerably 
broader.  The  multi-component  dispersion  relation  is  close  to  the  Vlasov 
result,  as  seen  in  Table  1. 

If  the  beam  is  propagating  into  initially  unionized  gas 
with  X  «  1,  Eq.  (14),  together  with  Eq.  (38)  or  Fig.  2,  indicates  very 
rapid  power-law  growth  of  the  instability.  If  on  the  other  hand,  the  gas 
is  pre-ionized  or  rapidly  ionized  at  the  beam  head,  the  limit  X  <<  1  should 
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correspond  to  propagation  in  a  fixed  channel.  This  limit  may  be  obtained 
by  writing 

=*  jc/;q  (39a) 

and  expanding  In  (r/r  )  in  Eq.  (14)  to  obtain  exponential  mode  growth  at  a 
rate  independent  of  X, 

A1(r,^)  =  A( r )exp[ -i/j)( ’  -  r  )/c].  (39b) 

Current  Neutralized  Case,  f  <  0,  Increasing  Values  of  1 

Figures  3-5  show  the  multicomponent  dispersion  relation  for  a  sequence 
of  increasing  values  of  X.  We  note  that  for  any  given  value  of  f  in  the 
range  -  1  <  f  <  0  (partial  current  neutralization),  Xui^  decreases 
as  X  increases.  Thus,  if  f  is  held  fixed,  the  growth  rate  jk  decreases 
with  increasing  1^  for  two  reasons:  m  *  X  1  *  1^  1  because  aQ(0  grows 
faster  for  large  1^,  and  additionally  because  for  large  1^  the  conductivity 
channel  tends  to  follow  the  perturbations  of  the  beam  current,  thus 
inhibiting  the  spatial  separation  of  1^  and  Ip.  The  latter  effect  is  most 
dramatic  at  3  «  0,  where 


-  f  (1  +  f)'1  -  X,  if  -  f  (1  +  f)'1-  X  >  o 
Xto^  ■  (40) 

0  ,  otherwise. 

For  moderate  values  of  f  and  large  values  of  X ,  X  'lax  ^  is  reduced  by 
about  30%  by  the  effect  of  perturbed  conductivity  over  the  range 
0  <  X  <  4.4  covered  in  Figs.  2-5.  In  the  limit  f  ♦  -  1,  however,  the 


stabilizing  effect  of  perturbed  conductivity  goes  away  and  x  \  ^  for 
fixed  f. 

The  spread  mass  dispersion  relation,  Eq.  (28),  differs  from  the 
multicomponent  results  mainly  in  that  the  unstable  spectrum  is  narrower  and 
the  peak  growth  rate  is  always  at  3  *  0.52  3  for  the  former,  as  seen  in 
Table  1.  The  Vlasov  unstable  spectral  width  is  close  to  the  multicomponent 
results. 

Ip  Parallel  to  1^,  f  >  0,  Increasing  Values  of 

The  effects  of  perturbed  conductivity  are  most  dramatic  if  f  >  0,  i.e. 
the  plasma  current  flows  parallel  to  the  beam  current  (as  may  occur  in  the 
back  half  of  a  triangular  beam  pulse  where  1^  is  falling,  or  in  the 
presence  of  a  strong  external  discharge).  Figure  2  shows  that 
when  X  +  0  the  hose  mode  is  damped  for  all  values  s>£  '{  if  f  >  C,7,  but 
Figs.  3-5  show  that  even  the  smallest  positive  \  leads  to  a  full  range  of 
unstable  modes,  and  that  if  X  »  1  the  external  discharge  does  not  even 
reduce  the  growth  rate  significantly.  (See  also  Table  I.)  In  this  limit 
the  plasma  current  tends  to  follow  the  distorted  beam,  rather  than 
constraining  the  beam  to  follow  a  straight  channel.  [Of  course  a  discharge 
may  also  exert  a  strong  stabilizing  effect  by  providing  a  pre-existing 
conductivity  at  the  beam  head;  mathematically  this  corresponds  to 
increasing  5  in  Eq.  (7).  This  effect  is  not  due  to  the  continuing 
presence  of  discharge  current  during  transport  of  the  beam. ] 

Limit  of  Large  X 

In  the  limit  of  very  large  X  and  fixed  f, 
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Eq.  (28)  reduces  to  the  dispersion  relation 


aX  -  iG(.i2/fl,  2),  (41) 

30 

identical  to  the  dispersion  relation  for  t^jjO)  in  a  fixed  channel  with  f  * 
0.  In  this  limit,  the  conductivity  channel  follows  the  hose  distortions  of 
the  beam  so  easily  that  Ip  does  not  separate  from  1^  and  presence  of 
equilibrium  plasma  current  has  no  destabilizing  effect  on  hose  stability. 
However  our  assumption  that  the  plasma  is  weakly  ionized  normally  will 
break  down  before  this  limit  is  reached. 

Beam  Injected  into  Neutral  Gas ,  No  Discharge  Current 

Of  particular  interest  is  the  scaling  of  hose  growth  with  1^ ,  for  the 
case  of  beam  propagation  in  an  initially  neutral  gas  with  no  externally 
driven  plasma  current.-  We  have  seen  that  the  hose  dispersion  relation 
depends  on  two  parameters:  X,  which  is  proportional  to  1^,  and  f,  which 
depends  in  general  on  the  radius  a  and  the  rise-time  Tr  of  1^(0,  as  well 
as  on  the  value  of  1^  in  the  beam  body.  However  a  simple  scaling  is 
possible  for  the  case  in  which  1^(5)  rises  quickly  to  its  plateau 

value  (t  / a  <  20  to  40),  and  avalanche  ionization  due  to  the  E  -spike  at 

f  v  6 

the  pinch  point  is  unimportant  (E  / p  <  100  keV/cm-atm  in  air,  i.e.  gas 

z  — 

density  not  too  low  and  beam  current  density  not  too  high).  Under  these 
conditions,  f($)  depends  only  on  1^;  after  reaching  a  maximum  at  the  pinch 

point,  it  levels  off  to  a  value  that  varies  slowly  with  ;  and  may  be 

estimated  as® 
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To  obtain  an  expression  in  closed  fora  for  the  scaling  of  the  peak  growth 

rate  with  1^,  we  use  (42)  and  (37)  in  the  spread  mass  dispersion  relation 

(28).  Noting  that  the  peak  growth  rate  Max  from  Eq.  (23)  always  occurs 

at  ft/ft„  9  0.52  where  G_  reaches  its  maximum  value  0.69,  we  find  that 
30  t 

Q 

Max  a  -  ■—  +  Re  [Q2  -  4Gf(l  +  f)'l/2  (43) 

with  Gr  *  0.69,  G^  *  1.01.  The  quantity  in  brackets  is  slowly  varying,  so 

that  Max  scales  essentially  as  f-^-.  Max  from  Eq.  (43)  is  plotted  as 

a  function  of  1^  in  Fig.  6,  using  the  multicomponent  model  and  Eq.  (42); 

the  result  differs  slightly  from  Eq.  (43).  We  see  that  Max  is  a 

monotonically  decreasing  function  that  approaches  a  limit 

Max  a.  *  0.69,  as  {  ^  *}.  (44) 

X  +  « 

Thus  higher  current  beams,  with  larger  return  current  fractions  |f|,  are 
seen  to  be  less  unstable. 

Equation  (42)  happens  to  be  the  condition  for  marginal  stability  of 
the  hose  mode  with  ft  *  0.  Since  (42)  is  only  a  rough  estimate  of  f,  and  f 
does  vary  somewhat  even  within  the  beam  body,  specific  cases  may  have  an 
unstable  mode  at  ft  •  0,  but  the  growth  rate  u^(O)  then  would  be  quite  small 
compared  to  Max  uj,  •  Instability  in  the  range  ft  <<  ft  is  driven  primarily 

1  (JO 

by  repulsion  between  1^  and  Ip  and  has  been  called  the  "self-hose",  but  as 
seen  in  Figs.  2-5  there  is  a  smooth  transition  to  the  much  stronger  "high- 
frequency  hose",  with  a  peak  growth  rate  in  the  range  0.6  <  ft/ft30  <  1. 

In  cases  where  avalanche  is  important  near  the  pinch  point,  or  where 
the  current  rise  time  is  long  If  I  can  be  larger  than  the  estimate  given 


by  Eq.  (42),  leading  to  a  Larger  hose  growth  rate.  Several  cases, 
calculated  with  the  multicomponent  simulation  code  VIPER,  (described  in  the 
next  section),  are  also  shown  in  Fig.  6.  For  these  cases,  the  relativistic 
factor  is  y  ■  100,  the  beam  current  pulse  shape  i3  V:)  -  *bo  tanh 
(;;/30ao),  the  gas  is  air  at  standard  density,  1^  is  4,  10,  40  or  ICO  kA, 
and  the  radius  of  the  beam  body  is  aQ  »  0.5  (1^/10  kA)  cm,  so  chat  the 
current  density  is  the  same  for  the  four  cases.  In  order  to  compare  with 
theory,  the  initial  perturbation  was  started  in  the  beam  body  at  ;  *  30aQ, 
and  recombination,  scattering  of  the  beam  and  ohmic  energy  loss  were  turned 
off.  The  hose  growth  rates  seen  in  Fig.  6  do  not  deviate  greatly  from  the 
scaling  of  Eq.  (43). 


3.  Instability  of  the  3eam  Tall 


As  '  Increases  and  the  plasma  electron  density  ng  becomes  large, 
recombination  can  become  competitive  with  impact  ionization,  and  thus 
“saturate"  the  plasma  conductivity  o.  This  effect  can  be  modeled  by 
extending  Eq.  (6)  to  the  form 


3_ 

3-; 


a0(r,;) 


<  Jbo<r>  - 


S  o 

TO 


(r). 


(45) 


For  simplicity,  we  neglect  temperature  dependendence  in  the  recombination 
coefficient  as  well  as  the  mobility.  For  a  weakly  ionized  gas  where 
electron-neutral  collisions  dominate  the  resistivity,  we  then  have 
3r  =  &T'  og/o0,  w^ere  is  a  constant  for  any  given  gas  species,  is 
the  gas  density,  and  pQ  is  standard  atmospheric  density;  a  typical  value 
for  air  is  Br*  7  x  10-15  sec/cm. 

The  solution  to  Eq.  (45)  is 


oo(r.O  -  <Jgat(r)  canh  (S^sat^’  (46a) 

where  the  saturated  value  is 

»8.t(r)  '  Jbo<r)/3rll/2-  (46b> 

We  see  from  Eqs.  (46)  that  a  saturates  first  on  axis,  where  <j  _  is 

o  sat 

largest;  saturation  can  occur  considerably  later  off  axis.  Of  course, 
saturation  may  not  occur  during  the  duration  of  the  beam  if  or  o  is  too 

small. 

We  consider  next  an  idealized  beam  tail,  where  oQ  is  saturated  at 


least  out  Co  several  beam  radii  a.  We  assume,  as  we  did  in  the  beam  body, 


that  I0  and  a  are  constant.  With  these  assumptions,  -  is  an  ignorable 
coordinate  in  the  beam  tail,  and  hose  eigenmodas  are  of  the  form  (22). 

The  perturbed  conductivity  7^  is  determined  from  the  linearized 
equation 


bl 


J0JL’ 


(47) 


i.e.  from  (22), 


23r3g(f)  “  ioi* 


(48) 


When  we  use  Eqs.  (48),  (46b)  and  (10a)  in  Eq.  (5),  the  linearized  Ampere's 
law  can  be  written  as 


dr  r  dr  2c 


iiroio  (r)a 

_2 - ] 

2XJbo(r)c  ' 


(49) 


analogous  to  Eq.  (15)  in  the  beam  body.  The  radial  profiles  of 

J  (r)  and  <7  (r)  are  essentially  identical  to  each  other,  but  are  broader 
po  o 

than  that  of  J,  (r),  so  the  radial  dependence  of  f(r)  =  J  (r)/J.  (r)  is 
bo  po  bo 

important. 

A 

As  in  Sec.  II,  J^( r )  for  use  in  Eq.  (49)  may  be  expressed  in  terms 

A 

of  A(r)  by  using  the  multi-component  model,  Eqs.  (30)  and  (31a),  or  the 
Vlasov  orbit  integral,  Eq.  (35).  For  the  multi-component  model,  the 


resulting  eigenvalue  equation  is 


where  r(r)  is  the  auxiliary  function  defined  in  Eq.  (33).  Equation  (50)  is 
analogous  to  Eq.  (32)  and  is  supplemented  by  Eq.  (31a)  and  the  boundary 
conditions  (34).  We  shall  use  Eq.  (50)  for  the  analysis  and  numerical 
examples  used  in  this  section. 

The  Vlasov  form  of  the  eigenvalue  equation  is  obtained  by  using  Eq. 
(35)  in  Eq.  (49).  This  integro-dif ferential  equation  may  be  solved 
numerically  by  our  code  VALIUM. 

In  regard  to  the  spread  mass  model  we  note  that  the  model  in  the  form 
of  Eqs.  (23)-(28)  can  only  be  used  if  J^Q(r)  and  oQ(r)  both  have  Bennett 
profiles  of  equal  width,  and  thus  is  not  applicable  to  the  present 
case.1^ 

Figure  7  shows  the  growth  spectrum  calculated  from  Eq.  (50)  for  four 
cases  spanning  a  wide  range  of  parameter  space:  X  ■  0.044,  f  -  -  0.042 
(solid  curve);  X  ■  0.44,  f  ■  -  0.30  (dashed  curve);  X  ■  4.4,  f  ■  -  0.81 
(dotted  curve);  and  X  ■  0.44,  f  ■  +  1.0  (dot-dashed  curve).  We  use  the 
notation  f  =  Ip/I^,  even  though  the  spatial  profiles  of  Jp(r)  and  J^(r)  are 
different.  The  first  three  cases  represent  roughly  beams  of  1,  10  and  100 
kA  Injected  Into  neutral  air,  to  the  accuracy  of  our  simple  conductivity 


formula  (45),  with  <  specified  by  (37)  and  our  return  current  scaling  given 
by  (42);  the  fourth  case  represents  a  10  kA  beam  "guided”  by  a  10  kA  net 
discharge  current. 

In  each  case  the  radial  profile  of  (r)  is  taken  to  be  a  3ennett 
profile  cut  off  at  three  3ennett  radii,  in  the  form 


For  this  profile,  the  Bennett  radius  as  is  slightly  smaller  than  the 

scale  radius  a  defined  in  Eq.  (2).  A  cut-off  is  necessary  because 

otherwise  the  slow  radial  fall-off  of  the  conductivity  profile  (46b)  would 

cause  essentially  all  of  the  plasma  current  to  flow  outside  the  beam.  This 

3  8 

is  one  of  many  mathematical  complications  ’  associated  with  the  slow  fall- 
off  of  the  Bennett  profile  at  large  radii;  the  3ennett  profile  is  a  proper 
representation  of  Jfeo(r)  only  in  the  central  core  of  the  beam.  A  cut-off 
is  also  consistent  with  the  slow  approach  to  conductivity  saturation  at 
large  radii.  In  cases  of  high  fractional  current  neutralization  the  growth 
rates  do  depend  on  the  choice  of  cut-off,  so  it  can  be  important  in 
modeling  experiments  to  determine  the  profile  of  J^0(r)  or  oQ(r)  in  the 
wings . 

It  is  striking  that  the  scaled  peak  growth  rate  x^  Max  is  nearly 

2 

the  same  (to  within  »  30%)  in  all  cases.  But  =  toq( 0)  a  /2c  increases 
rapidly  with  1^;  if,  for  example,  the  beam  current  density  J^(0)  on  axis  is 
held  constant  as  1^  is  varied,  then  the  recombination  saturated  value 
of  oo(0)  is  constant  while  a  *  Ib  and  x^  «  1^.  Thus  Fig.  7  shows  that  the 
actual  growth  rate  in  the  beam  tail  Max  u>^  decreases  as  even  though 

the  return  current  fraction  increases  rapidly  with  Iv>. 


is  due  to  three 


The  strong  decrease  of  Max  with  increasing 

separate  effects:  (i)  The  conductivity  channel  relaxes  toward  the 

perturbed  beam,  thus  reducing  the  magnetic  repulsion  between  the  plasma 

return  current  and  the  beam.  This  effect  appears  in  Eq.  (47)  and  thence  in 

the  f(r)  term  in  the  first  factor  of  the  right  hand  side  of  Eq.  (49). 

(ii)  The  unperturbed  conductivity  channel  .7  (r)  is  broader  than  the  beam 

profile  'Jbo(r)‘  as  seen  *-n  ^9*  (46b).  Thus  much  of  the  plasma  current 

flows  outside  the  beam  and  exert3  no  force  on  the  beam.  Mathematically, 

the  destabilizing  factor  [1  +  f(0)]”^  on  the  right  hand  side  of  Eq.  (49)  is 

smaller  (for  Ip/1^  <  0)  than  (1  +  Ip/I^)-^.  (iii)  The  growth  rates 

naturally  scale  inversely  with  the  dipole  diffusion  time,  some  radial 

2  2 

average  of  the  quantity  iraQ(r)a  /2c  that  appears  on  the  left  hand  side  of 
Eq.  (49).  When  Jtj0(r)  an£*  aQ(r)  **ave  identical  Bennett  profiles,  the 
appropriate  average  is  t^,  defined  in  Eq.  (4),  but  when  the  profile 
of  o0(r)  is  broader  than  that  of  J|,0(r),  the  true  dipole  diffusion 

^  a  1  1  m 

time  x^  is  longer  than  x^5  this  effect  reduces  Max  for  all  cases,  ~  * 

even  if  1 1  I  «  I,  . 

P  a 

The  numerical  results  of  each  of  these  effects  is  illustrated  in  Table 
2.  For  each  of  the  four  cases  calculated  in  Fig.  7,  x^Max  as  calculated 
for  the  correct  equilibrium  (first  column)  is  compared  with  the  value 
obtained  if  the  conductivity  channel  is  not  allowed  to  follow  the  beam 
perturbations  (second  column),  or  if,  in  addition,  the  channel 
shape  <jQ(r)  is  set  equal  to  beam  shape  d^Q(r).  We  note  that  the 
stabilizing  effects  are  present  in  all  cases,  but  are  very  strong  for 
highly  current-neutralized  cases. 


■  ->  ”> 
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Table  2 

Comparison  of  the  non-dimensionalized  peak,  growth  rate  Max  ^  in  the  beam 
tail,  as  calculated  from  the  correct  equilibrium  [rQ(r)  *  i^“(r)].  with 
perturbed  conductivity  j  included,  with  j  neglected,  and  from  a  prescribed 
fixed  conductivity  channel  with  a^Cr)  <*  J^Cr).  Four  cases  are  considered: 
the  first  three  are  approximately  correct  values  of  f  for  a  beam  with  the 
given  y  injected  into  neutral  gas,  and  the  fourth  case  is  for  a  plasma  current 
externally  constrained  to  be  equal  to  1^. 


Maximum  Growth  Rate  Max 


1/2 


°o  *  Jbo  '  *  Jbo 


i/: 


a  self-consistent  a  *  0 


7o  Jbo 
a  ■  0 


0.044 


-0.04 


0.49 


0.49 


0.63 


0.44 


-0.3 


0.53 


0.58 


1.20 


4.4 


-0.8 


0.40 


0.81 


7.0 


0.44 


0.40 


0.25 
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Multi-Component  Simulations 


A.  VIPER  Model  Description 

The  VIPER  multi-component  simulation  model  allows  one  to  relax 
many  of  the  simplifying  assumptions  of  the  previous  sections.  In 
particular,  this  code  can  treat  the  beam  head  region  where  the  assumptions 
of  constant  I^(^),  I^(;),  a(?)»  and  aQ(c)  are  violated,  and  where  the 
magnetostatic  approximation  to  the  electromagnetic  fields,  Eq.  (5),  is  not 
valid.  Because  the  conductivity  in  this  region  is  low,  considerable  hose 
growth  can  occur  before  the  beam  body  is  reached.  In  addition,  simulations 
indicate  that  the  dominant  frequency  2  observed  in  the  beam  body  and  tail 
is  strongly  influenced  by  conditions  in  the  beam  head.  The  sharp 
separation  of  the  beam  into  a  body  and  a  tail,  as  described  in  the 
dispersion  relation  models,  is  also  somewhat  arbitrary  since  conductivity 
saturation  due  to  recombination  occurs  first  on  axis  and  much  later  in  the 
radial  wings  of  the  beam.  Thus,  there  is  a  long  transition  region 
in  t,  between  the  body  and  the  tail. 

The  VIPER  model  treats  the  evolution  of  self-consistent  fields, 
conductivity,  and  beam  dynamics  for  both  monopole  (m-0)  and  dipole  (m»l) 
quantities.  A  description  of  the  model  equations  is  given  in  the 
Appendix.  The  coordinate  transformation  to  r,  and  z  discussed  in  the 
introduction  is  used  throughout,  and  all  quantities  are  radially 
resolved.  Monopole  and  dipole  fields  are  treated  using  the 
ultrarelativistic  equations  derived  by  Lee  for  arbitrary  conductivity. 

In  addition  to  the  beam-impact  ionization  term  in  (6),  the  conductivity 
equations  include  a  simple  E/p  avalanche  model  and  recombination. 
Axisyoaetric  beam  dynamics  are  represented  by  an  envelope  equation  which 
assumes  that  the  beam  expands  or  contracts  self-similarly ,  maintaining  a 
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Bennett  profile.  Zmittance  changes  due  to  multiple  scattering  and 
anharmonic  damping  influence  the  envelope  radius,  as  does  beam  energy  loss 
due  to  the  inductive  electric  field  EZQ.  Dipole  beam  dynamics  are  treated 
using  the  linearized  multi-componer.c  model  described  in  Ref.  3.  Each  axia 
segment,  propagating  at  a  fixed  distance  r  from  the  beam  head,  is 
partitioned  into  several  hundred  components  whose  densities  are  chosen  so 
as  to  reproduce  the  desired  Bennett  profile.  Since  the  individual 
components  have  different  edge  radii  and  since  each  is  allowed  to  oscillat 
independently  in  the  transverse  plane,  the  dynamics  are  radially 
resolved.  The  principal  quantity  of  interest  is  the  average  hose 
displacement  Y(^,z)  of  these  components,  weighted  by  the  individual 
component  densities.  Because  the  model  is  linearized,  Y  is  treated  as  an 
infinitesimal  (compared  to  the  beam  radius),  and  is  normalized  to  the 
initial  perturbed  displacement. 


B.  VIPER  Simulation  :>f  a  10  kA  Beam 


In  this  section,  we  present  an  example  of  a  single  VIPER 
simulation  of  a  50  MeV,  10  kA  beam  propagating  in  full  density  air. 
Parameters  for  the  run  are  listed  in  Table  3.  All  distances  are  expressed 
in  terms  of  the  nominal  beam  radius  aQ  ■  0.5  cm. 

The  beam  radius  a,  on-axis  conductivity  ,  and  "effective  current 

fraction"  1  +  fg  are  plotted  as  functions  of  distance  ;  from  the  beam  head 

in  Fig.  3.  Here  1  +  f  is  the  pinch  force,  averaged  over  the  bean  profile 

J^Q(v),  and  normalized  to  the  value  it  would  have  if  there  were  no  plasma 

current  and  the  self-fields  were  purely  magnetostatic.  Thus  fg  =  Ip/Iv, 

the  beam  body,  where  dpo(r)/Jk0(r)  is  approximately  constant,  but  1  + 

fe  +  0  in  the  beam  head,  where  space  charge  forces  largely  cancel  the 

magnetic  pinch  forces,  and  in  the  beam  tail  the  difference  between  the 

profiles  of  Jpo(r)  and  J^^r)  *s  *-ncluded  in  the  definition  of  fg.  The 

beam  radius  a(^),  plotted  at  five  values  of  the  propagation  distance  z/aQ 

7  8 

between  0  and  2400,  shows  the  processes  of  nose  erosion  *  and  Nordsieck 
expansion^  (Fig.  8a).  The  conductivity  r~0),  seen  in  Fig.  8b,  rises 

linearly  with  5  after  1^  reaches  its  plateau  value  (?/aQ  >  50),  and  until 
recombination  saturation  sets  in.  As  a  function  of  z,  aQ  decreases  slowly 
due  to  the  effect  of  Nordsieck  expansion.  The  effective  current  fraction 
(Fig.  8c)  varies  slowly  with  both  5  and  z  for  >  50.  Thus,  this  beam 

has  a  substantial  region  where  a,  Ip/I^,  and  daQ/d^  are  approximately 
constant  in  5  and  z,  and  the  dispersion  relation  analysis  for  the  beam  body 
can  be  applied. 

Figure  9  plots  the  average  hose  displacement  Y(z)  for  ;/a  *  124 

and  149.  Because  the  instability  is  convective  in  the  beam  frame,  Y(z)  at 
a  fixed  location  t,  initially  grows,  but  eventually  saturates  and  decays. 
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Table  3 


Parameters  for  VIPER  Simulation  Example 


Parameter 


Symbol 


Value 


Maximum  3eam  Current 

Initial  relativistic  factor 

Nominal  radius 

Gas  pressure 

Rise  length 

Pulse  length 

Start  of  perturbation 

Initial  perturbation  form 

Current  profile 


P 


10  k.A 

10C 


0.  5  cm 
760  torr 


*r 

r 

’’max 

V«> 


300ao 

20ao 

sin  [0. 157(5-ro)/ao J  for  20  <  */a  <  40 

!b0  canh  (?/?r) 
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The  frequency  ?.  of  the  hose  oscillations  seen  in  Fig.  9  varies  slowly 
with  r  and  2,  so  a  single  mode  approximation  in  the  dispersion  relation 


(28)  or  (32)  can  be  used.  Taking  local  beam  parameters  at  ;/ a  =  119  and 

o 

z/aQ  *  600,  we  have  -l-^a  ■  0.057,  and  o/f!  =  0.31.  Note  that  this  value 

is  significantly  below  the  value  ?.!?.  =  0.62  expected  from  Eq.  (32)  and 

Table  1  for  peak  hose  growth. 

For  a  convective  instability,  the  maximum  ("saturated")  value 

of  Y  reached  at  a  given  ”,  denoted  Y  (;),  should  increase  with  r  at  a 

S3 1 

power-law  growth  rate  in  the  beam  body,  as  prescribed  in  Eq.  (14)  or  at 
an  exponential  growth  rate  ^  in  the  beam  tail,  as  prescribed  in  Eq. 

(22).  For  example,  |Y  |  (^)  is  1900  for  r  =  124,  but  is  3000  for  ;  * 

SaL 

149,  as  seen  in  Fig.  9.  The  plot  of  Yga  (^),  Fig.  10,  shows  power  law 

growth  as  expected.  A  least  squares  fit  gives  u)^exp  =■  2.28  for  the 

"computer  experiment"  growth  rate.  The  multicomponent  dispersion  relation 

for  1  +  f  *  0.57,  X  ■  0.44,  and  ?-/2^Q  *  0.31  gives  *  2.40. 

(Variations  if  f  and  Cl/ Cl  within  the  beam  body  introduce  an  uncertainty 

So 

of  +  0.25  in  the  theoretical  prediction  for  w^.)  The  observed  growth  rate 

thus  agrees  well  with  the  dispersion  relation  predictions  for  the  observed 

value  of  Cl,  but  is  substantially  below  the  value  =  3.8  predicted  for 

the  fastest  growing  mode,  where  ft/fl  *  0.62. 

This  reduction  in  hose  growth  rate  from  the  theoretical  maximum  is 

seen  in  all  VIPER  simulations  in  which  the  perturbation  begins  in  the 

expanded  beam  head.  This  region  of  the  beam  has  a  lower  Q  than  in  the 
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body  and  a  lower  dipole  decay  length.  Low  frequency  modes  grow  rapidly  in 
the  beam  head  and  persist  into  the  beam  body.  Higher  frequency  modes  are 
stable  in  the  beam  nose  and  usually  do  not  catch  up  with  the  lower 
frequency  modes  in  the  beam  body.  Thus,  the  observed  Cl/ Cl^Q  in  the  beam 
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he  amount  of  hose 


body  is  substantially  less  than  0.62.  In  general,  t 
growth  depends  strongly  on  the  extent  to  which  -/arias  throughout  the 

beam. 

In  the  example  used  here,  the  simplifying  assumptions  of  the 

analytic  theory  of  Sec.  II  are  very  well  satisfied  within  the  beam  body. 

Even  in  this  case,  we  have  seen  that  the  transition  from  beam  head  to  body 

plays  an  important  role  in  choosing  the  dominant  value  of  2.  In  many  other 

cases,  the  assumptions  may  not  be  well  satisfied,  so  that  the  analytic 

theory  is  only  qualitatively  applicable.  For  example,  the  I^(;)  pulse 

shape  mav  rise  slowly  or  may  not  have  a  flat-too.  I  (-)  nav  fall  off 

P 

significantly  within  the  beam  body,  and  therefore  the  beam  radius  may 
decrease  as  the  pinch  becomes  stronger,  particularly  for  high-current  cases 
where  -  1  <  I^/I^  <  -  -j.  In  such  cases,  code  studies  are  essential  to 
predict  hose  growth  accurately.  A  variety  of  such  studies,  over  a  wide 
parameter  range,  will  be  presented  in  future  reports. 


5. 


Conclusions 


We  have  derived  dispersion  relations  that  include  conductivity 
evolution  self-consistently  (but  with  the  important  assumption  of  constant 
electron  mobility)  in  both  the  equilibrium  and  the  linearized  hose 
perturbations  of  a  propagating  beam.  We  have  used  the  dispersion  relations 
to  explore  a  wide  range  of  beam  and  gas  parameters.  Particularly  for  cases 
in  which  there  is  a  high  degree  of  overall  current  neutralization,  the  hose 
growth  rate  turns  out  to  be  much  less  than  might  have  been  expected  from 
previous  work.  We  have  also  used  the  multicomponent  code  VIPER  to 
demonstrate  some  additional  effects  that  can  modify  the  interpretation  and 
use  of  the  dispersion  relations.  As  an  example,  hose  growth  in  the  beam 
body  is  usually  dominated  by  low  frequency  modes  that  are  excited  at  and 
grow  rapidly  near  the  pinch  point  at  the  front  of  the  beam,  rather  than  by 
the  modes  that  are  most  unstable  in  the  body. 
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Appendix 


In  this  appendix  we  summarise  the  model  equations  used  in  the  VIPER 
hose  dynamics  code  discussed  in  Section  4.  The  model  is  based  on  the 
linearized  multicomponent  formulism  described  in  Ref.  3.  Ve  represent  the 
unperturbed  beam  current  J^Cr)  ns  a  superposition  of  components  with 
density  profile 


(1  -  r2/R2)  G(R  -  r), 


(A-l) 


where 


0(R  -  r)  s 


1,  if  J  )  r, 
0,  if  R  <  r. 


The  edge  radius  R  is  used  to  index  the  components.  Thus  we  can  write 


2,. 2, 


Jfeo(r)  -  /  dR  <1  -  r  /R  ;  F(R)  9(R  -  r), 


(A-2) 


where  F(R)  is  the  amplitude  of  component  R.  F(R)  can  be  calculated  by 
differentiating  Eq.  (A-2)  twice. ^  For  a  beam  with  3ennett  profile 


Jbo(r’  z) 


iraZ(5,z)  [1  +  r2/a2(5, z)  ] 2 


(A-3) 


we  find 


?(R)  - 


121.  R' 
0 


*a6(l  +  R2/a2)4 


(A-4) 


.13 


The  perturbed  current  Jvi(r)e  is  calculated  by  assuming  that  each 
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component  undergoes  a  rigid  linearized  displacement  Y(R,s,z)  along  the  y- 
axis  (perpendicular  to  the  direction  of  propagation).  The  formal  result  is 


Jbo(r)  +  Jbl(r)sin9 


{  dR  9  [R  -  | r  -  Y(R) I ”]{ 1  - 


ir  -  Y(R) 


F(R), 


(A- 5) 


To  take  account  of  the  fact  that  a  slice  may  expand  or  contract,  we 
generalize  the  formalism  slightly  by  defining  p  =  R/a,  and  use  p  as  the 
component  index,  i.e.  each  component  expands  at  the  same  rate  as  a  and  the 
beam  expansion  is  thus  self-similar.  Linearizing  in  Y(R,;,z)  and  expanding 
the  theta  function  gives 


Jbl(r) 


r 

r  / 


F(p?  Y(p) 

2 

P 


24I.r  . 

b  f  »  oY(p) 

- -  j  dp  - j-j 

ia  r/a  (1  +  p“) 


(A-6) 


for  a  beam  of  Bennett  profile.  In  practice,  the  VIPER  model  includes  only 
components  with  p  <  Pmair"  3;  as  a  result,  the  unperturbed  current  J^0(r), 
which  is  calculated  by  evaluating  (A-2)  for  p  <  pmax»  falls  off  more 
rapidly  than  a  true  Bennett  distribution  for  r/a  >  3. 

In  VIPER  we  normally  use  a  conductivity  evolution  equation 


do 

d? 


*  Jb 


(A-7) 


which  includes  avalanche  as  well  as  the  beam-impact  ionization  and 
recombination  included  in  Eq.  (45).  (We  have  used  more  elaborate 
temperature-dependent  conductivity  models  on  occasion.)  The  avalanche 
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coefficient  is  a  function  of  the  ratio  of  the  electric  field  ZQ 
to  the  gas  density  o  ,  given  as  Eqs.  (16)  and  (17)  of  Ref.  13.  In  cgs 

O 

units,  the  recombination  coefficient  used  in  VIPER  for  air  is 

3  •  7.1  x  10  15(o  !q  )  sec/cm,  where  o  is  the  air  iensitv  at  standard 
r  g  o  o 

temperature  and  pressure,  and  the  impact-ionization  coefficient  is  <  •  3.6 
x  10“^cm/statamp-sec.  An  appropriate  expansion  of  (A-7)  gives  da^/d^. 

The  electromagnetic  fields  are  calculated  using  Lee's  "Empulse" 
approximation*-  to  Maxwell's  equations,  which  is  appropriate  for 
ultrarelativistic  paraxial  beams.  In  the  Lorentz  gauge,  these  field 
equations  are 


c_  l  3_  3_  f  *  , 

4ir  r  3r  r  3r  o  o  • 


bo 


+  a 

o 


* 

3  A 

o 


(A-8) 


c_  1_  3_ 
4ir  r  3r 


r 


2  * 
3  A 

o 

3r3c 


1  3 _ 

r  3  r 


ra 


o  3r 


(A-9) 


for  the  a  ■  0  potentials.  We  have  defined 


A  =  A  —  4  . 
0  zo  r0 


(A-10) 


The  corresponding  a  •  1  equations,  linearized  in  perturbed  quantities,  are 


£_  1_  I  rrA  ’ 

4ir  3r  r  3r  ^  1 


.  ,  „  3  A, 

c  3  1  3  _  _ 1_ 

4ir  3r  r  3r  3? 


♦  xJ  "  - 

3  A. 

3  A 

0 

Jbl  ao  3? 

+  «T1 

3?  ’ 

1  3 

r  3r  r 

3<f>. 

^o—  +al 

34 

_£)  - 

3  r; 

Vi 

*• 

(A-ll) 


(A-12) 


where 


Within  Lee's  approximation  the  electric  field  components  are  given  by 


3A 


E.  -  -  ?,p, 

-1  i 


:a-u) 


(A-15) 


and  the  Lorentz  force  on  the  beam  electrons  is  given  by 


FLz  "  "  eV 


F  -  e7  A  . 
-Li  1 


(A-16) 


(A-17) 


The  boundary  conditions  used  are  perfectly  conducting  conditions  on  a 

cylinder  at  large  radius,  as  well  as  the  conditions  that  all  fields  vanish 

ahead  of  the  beam*  These  field  equations  are  valid  in  low  conductivity 

regions  as  well  as  the  high  conductivity  beam  body. 

Axisymmetric  beam  dynamics  are  treated  using  an  envelope  model  of  a 

type  widely  used  in  beam  physics  problems. We  define 
2-22  -2 

n  "a  /a  ,  where  a  is  the  mean  square  radius  of  the  beam  (n  depends 
weakly  on  the  cut-off  rMX  applied  to  the  Bennett  profile^’®),  and 

2  7 

determine  the  constant  n  by  integrating  r‘J^0(r)  out  to  the  beam  edge 
rmax*  resulting  equations  describing  the  m  ■  0  dynamics  are: 


2  2  2  2  23  da 

d  a  _  e  -yaU-yna  dz  dz 

77 - — nrr^ - * 

dz  yap 


(A-18) 


yy?*.*  ?  ?*y  7  '■■  *.•  7  ?  *.-  i.^i.-v,.p|,^.'’'.-'! 


dc 


dz 


2a  ,na^y  ^;L’(d“a/ dz“) 
a 


7  7 

y'a'U 


+ 


( A- 19) 


dy  _  e  „ 

dz  ^  zo" 
me 


(A- 20) 


2 

Here  e  is  the  squared  beam  emittance,  U  is  the  radially  averaged  pinch 

*■>  9 

force,  is  a  sausage  oscillation  damping  coef f icient““  (taken  to  be  0.7) 
Ezq  is  the  axial  electric  field  at  r*0,  and  z^^  is  the  multiple  scattering 
rate. 


s  ■  (2h  c/e"  X  ) 
sc  r 


(A-21) 


where  h  is  Planck's  constant  and 


X_-  2.7  x  104(p /p  )ca 
r  o  § 


(A-22) 


is  the  radiation  length  in  air.  The  pinch  force  term  is 

0 - ^ 

■'Vc 

Finally,  m*l  dynamics  are  described  by  equations  for  the  hose  displacement 
Y(p)  of  components  at  a  fixed  value  of  ?.  The  radially  averaged  y 
component  of  force  on  the  beam  component  denoted  by  p  is  found  to  be 


»  ,  dA*(r) 

U dr  r  Jbo^>  -V' 


(A-23) 


4  F(p)  aAp  ap  -  dA  . 

- 3-/  dr  r2[Y(,)-^f+Al 

(ap)4  0  dr  1 


(A-24) 


where  Ap  is  the  spacing  between  components  (an  infinitesimal  quantity  in 
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analytic  expressions) .  A  Newton's  law  equation  of  motion  for 

Y(p,  S,  z)  can  be  obtained  by  dividing  out  the  component  mass  density, 

giving 


9z  (ap)  me  0 


ao  -  *  dA 

-  f  dr  r2  [\  +  Y  -~ 

4  -  J  '  -  1  dr  • 


(A-25) 


A  more  detailed  derivation  of  the  above  equations  and  a  description  of 
numerical  methods  is  found  in  Refs.  9. 
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Figure  1 

The  complex  function  G  defined  in  Eq.  (26),  as  a  function  of  x  = 
(upper  scale)  or  of  n/3.  (lower  acale). 
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Figure  2 

Growth  rates  in  the  beam  body  for  X  ■  0  and  several  values  of  f. 
Calculations  are  based  on  the  multi-component  model,  Eqs.  < 3 2) — ( 34) . 

Note  that  the  ordinate  is  (1  +  f)\oj^,  a  form  which  conveniently  compress 
the  range  of  variation. 


Figure  3 

Growth  rates  in  the  beam  body  for  \  ■  0.044  and  several  values  of  f 
Calculations  are  based  on  the  multi-component  model,  Eqs.  (32) — ( 34 ) 
Note  that  the  ordinate  is  (1  +  f),\jj. . 
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Figure  4 


Growth  rates  in  the  beam  body  for  \  *  0.44  and  several  values  of  f. 
Calculations  are  based  on  the  multi-component  model,  Eqs.  (32) — ( 34) . 
Note  that  the  ordinate  is  (1  +  f)Xu^. 


Growth  rates  in  the  beam  body  for  \  ■  4.4  and  several  values  of  f. 
Calculations  are  based  on  the  multi-component  model,  Eqs.  (32)— (34) 
Note  that  the  ordinate  is  (l  +  f).\y, . 
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Figure  6 

Maximum  growth  rate  in  the  beam  body,  for  a  beam  injected  into  neutral 
gas,  as  a  function  of  \  (lower  scale)  or  Ib  (if  the  gas  is  air,  upper 
scale).  The  curve  is  the  theoretical  result  from  Sqs.  (42)  and  (43), 
while  the  dots  are  the  results  of  four  VIPER  runs. 
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multicomponent  model,  Kq.  (50). 


Beam  equilibrium  properties  from  the  VIPER  rur  described  in  Table  3. 
Beam  radius  a/aQ,  (b)  conductivity  <jo(r  »  0),  and  (c)  pinch  strength 
Ie/I|)  s  1  +  f  ,  each  at  a  sequence  of  different  propagation 


ranges  z 
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13.  The  particular  expression  (4)  applies  for  the  case  of  identical 


's! 
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